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Summary 


In  the  first  part  of  this  research  a  general  geometrically-exact  Lagrangian  me¬ 
chanics  formulation  for  the  aerodastic  analysis  of  a  maneuvering  aircraft  has 
been  presented#  see  Morino  and  Baillieul  (1988)  for  the  details.  The  motion  of 
the  aircraft  is  expressed  in  terms  of  the  location  of  the  origin  of  a  body  frame  of 
reference,  the  rigid-body  rotation  of  the  body  frame  of  reference,  and  a  defor- 
mation.jThe  orthogonal  rotation  matrix  describing  the  orientation  of  the  body 
frame  with  respect  to  a  chosen  inertial  frame  of  reference  is  used,  aa  Lagrangian 
coordinate  for  the  rotational  .degrees  of  fi^edoinV^jrhe  Lagrangian  equations  of 
motion  for  the  corresponding  degrees  of  freedom  have  been  obtained. 

In  the  second  paxt  of  this  research^  formulation  has  been  specialized  to  the  case 
of  a  fluttering  buckled  plate  on  an  aircraft  undergoing  a  pitching  maneuver#see<- 
Sipcic  and-Morino  (1989)  for  detailtf.^  Assuming  that  the  maneuvering  of  the 
aircraft  is  prescribed,  the  Lagrange  equations  of  motion  for  the  elastic  degree  of 
freedom  has  been  derived.  This  equation  is  then  used  to  study  the  response  of 
the  panel  of  an  aircraft  engaged  in  a  pull-up  maneuver.  The  large-amplitude 
responses  are  investigated  by  using  the  digital  computer.*^  Long-time  histo¬ 
ries,  three-dimensional  view  of  orbits,  phase  planes,  and  power  spectra  of  the 
response  are  presented.  The  numerical  simulation  for  periodic  and  chaotic  re- J 
sponse  are  conducted  in  order  toCanalyze  the  influence  of  the  maneuvering  to  the 
dynamic  behaviour  of  the  panel)“>As  the  maneuvering  (load  factor)  increases, 
system  exhibits  complicated  dynamic  behavior  including  period-multiplying 
and  demultiplying  bifurcations  and  chaos.\ 
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1  The  Equation  of  Motion 

Let  us  summarize  the  formulation  and  the  results  of  the  first  part  of  this 
research  presented  in  Morino  and  Baillieul  (1988);  details  of  the  formulation 
are  in  the  original  report.  As  mentioned  before  the  motion  of  the  aircraft  is 
expressed  in  terms  of  the  motion  of  the  origin  of  a  body  frame  of  reference,  the 
rigid-body  rotation  of  the  body  frame  of  reference,  and  a  deformation,  which  in 
turn  is  expressed  (in  the  body  frame  of  reference)  in  terms  of  specified  functions 
^r(^“)i  as 

«((“.<)  =  E“-w.(n  (1) 

r 

Introducing  an  unconventional  approach  (t.e.,  using  the  rotation  matrix  as  a 
Lagrangian  coordinate),  and  the  following  the  conventional  procedure  for  La- 
grangian  mechanics  (t.e.,  calculating  the  kinetic  energy,  the  potential  energy, 
and  the  work  done  by  the  external  forces,  and  substituting  these  in  the  La¬ 
grangian  equations  of  motion)  one  obtains:  for  the  translational  degrees  of 
freedom, 

mOtXo  +  Dj  •  =  m(vo  +  Ovo)  +  hso  -I-  n*to 

+  ^  (iti,  +  2n«r  -I-  ftUr  +  S,.  =  fo  =  mR''  g  +  Cfl  (2) 

r 

for  the  rotational  degrees  of  freedom, 

Dtho  +  s  X  DeVo  =  mx  +  uirtti.  +  •  x  gB  (3) 

and  finally  for  the  elastic  degrees  of  freedom, 

•r  +  bj  «  -H  -  2w'' 

8£ 

+  ^  =  (4) 

For  the  definition  of  the  coefficients  in  terms  of  the  as  well  as  the  details 
of  the  derivation  the  reader  is  referred  to  the  original  report. 

The  formulation  includes  all  the  geometric  non-linearities,  but  has  been  fo¬ 
cussed  primarily  on  the  dynamic  ones,  which  arise  from  the  coupling  between 
the  ripd-body  rotation  of  the  frame  of  reference  and  the  other  degrees  of  free¬ 
dom.  The  structural  ones  arise  from  the  non-linear  terms  in  the  expression  of 
the  strain  tensor. 
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Finally,  the  simplification  resulting  from  the  use  of  the  natural  modes  of  vi¬ 
bration  for  the  functions  4,  has  been  indicated.  The  relationship  between 
the  finite-element  and  vibration-mode  formulation  has  been  examined.  It  was 
shown  that  finding  the  approximate  modes  of  vibration  from  a  finite-element 
formulation  and  using  the  corresponding  Lagrange  equation  of  motion  is  fully 
equivalent  to  transforming  the  system  of  finite-element  equations  in  such  a 
way  that  the  linear  portion  of  the  system  is  diagonal.  If  the  natural  modes 
of  vibration  is  less  than  the  number  of  unknowns  in  the  original  finite-element 
equations,  this  corresponds  to  reducing  the  number  of  unknowns  (and  the  cor¬ 
responding  equations)  in  the  transformed  system. 

2  Dynamic  of  Fluttering  Plate  in  Maneuvering  Aircraft 

In  order  to  study  the  effect  of  maneuvering  on  the  elastic  response  of  the 
fluttering  plate  in  the  maneuvering  airplane  we  assume  the  rigid-body  motion 
of  the  frame  of  reference  to  be  prescribed.  Hence  Eqs.  2,  3  are  not  used, 
and  in  Eq.  4  w  is  prescribed  and  therefore  the  inertial  terms  are  linear  in 
Ur  (with  time  dependent  coefficients).  A  constant  velocity  pull-up  maneuver 
has  been  assumed.  Following  the  classical  panel— flutter  approximation  and  the 
conventional  procedure  of  Lagrangian  mechanics,  the  equations  of  motion  are 
pven  by 

u  -f-  fu  -t-  Gu  +  tid  +  fa>v  =s  p  -f  f(tt)  7u»*u  (5) 

For  the  definition  of  the  coefficients,  as  wdd  as  the  details  of  the  derivation  the 
reader  is  referred  to  the  original  report. 

Equations  5  are  a  set  of  ordinary,  nonlinear  differential  equations  in  time.  The 
cubic  type  of  non-linearities  are  of  a  geometric  origin  and  are  associated  with 
the  occurrence  of  tensile  stresses  in  the  middle  surface.  The  coupling  between 
the  rigid-body  rotation  of  the  frame  of  reference  and  the  elastic  degrees  of 
freedom  is  represented  by  the  fourth  term.  Note  that  by  assuming  w  =  0  Eqs. 
5  properly  reduce  to  a  classical  panel  flutter  equations. 

S  Numerical  Experimentation 

The  extensive  numerical  simulation  for  periodic  and  chaotic  response  are  con¬ 
ducted  in  order  to  analyze  a  behaviour  of  the  fluttering  panel  of  the  maneuver¬ 
ing  airplane,  for  details  see  oripnal  report  by  Sipdc  and  Morino  (1989).  Several 
important  conclusions  derive  from  the  preceding  work,  see  Fig.  1  and  sequence 
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of  Figs.  2(a)  thru  (r)  as  the  example.  The  new  limit  cycle  has  been  observed 
in  the  non-maneuvering  case.  It  was  shown  that  chaos  could  occur  in  a  ma¬ 
neuvering  case  for  system  parameters  in  the  actual  flight  range.  The  presence 
of  a  load  factor  can  transform  the  response  from  the  fixed  point  into  a  simple 
periodic  or  even  chaotic  state.  It  can  also  suppress  the  periodic  character  of  the 
motion,  transforming  the  response  into  a  fixed  point  (compare  responses  Fig. 
2(a)  and  (r)).  Regions  of  the  chaotic  motion  with  an  island  of  non-sinusoidal 
periodic  motion  has  been  observed,  see  Figs.  2(1)  thru  2(q).  Relatively  sudden 
jumps  between  these  two  regimes  were  found.  For  some  values  of  the  load  fac¬ 
tor,  a  cascade  of  period-multiplying  bifurcations  begins,  culminating  in  chaos 
Figs.  2(d)  thru  2(h),  followed  by  a  series  of  period-demultiplying  bifurcations 
as  the  system  returns  to  a  starting  limit  cycle  response  Figs.  2(h)  thru  2(k). 
The  region  of  intermittent  and  transient  chaos  has  been  observed  as  well. 

4  Concluding  Remarks 

A  general  geometricaUy-exact  Lagrangian  mechanics  formulation  for  the  aeroe- 
lastic  analysis  of  a  maneuvering  aircraft  has  been  presented.  The  motion  of 
the  aircraft  is  expressed  in  terms  of  the  location  of  the  origin  of  a  body  frame 
of  reference,  the  rigid-body  rotation  of  the  body  frame  of  reference,  and  a 
deformation,  given  as  a  linear  combination  of  prescribed  functions,  The 
Lagrangian  equations  of  motion  for  the  corresponding  degrees  of  freedom  have 
been  obtained.  The  advantages  of  identifying  4r  natural  modes  of 

vibration  (instead  of  finite-element  test  functions)  have  been  discussed. 

In  the  second  part  of  this  research  a  general  geometrically-exact  Lagrangian 
mechanics  formulation  for  the  aeroelastic  analysis  of  a  maneuvering  aircraft 
has  been  specialized  to  the  case  of  a  fluttering  plate  undergoing  a  pitching 
maneuver.  The  formulation  includes  geometric  non-linearities  associated  with 
the  occurrence  of  tensile  stresses  in  the  middle  surface,  as  well  as  coupling  be¬ 
tween  the  rigid-body  rotation  of  the  frame  of  reference  and  the  other  degrees 
of  freedom.  The  general  response  of  the  system  were  simulated  on  the  digi¬ 
tal  computer  by  using  a  fourth  order  Runge-Kutta  algorithm  with  adaptive 
step  size  control.  Long-time  histories,  three-dimensional  view  of  orbits,  phase 
planes,  and  power  spectra  have  been  used  to  characterize  the  response. 

The  results  indicate  that  the  study  of  this  deterministic  system  is  important 
from  the  practical  and  theoretical  viewpoints.  First,  the  fact  that  maneuvering 
can  change  the  character  of  the  panel  response  is  of  practical  interest  to  the 
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aeroelasticians  as  it  affects,  for  instance,  fatigue  analysis.  Secondly,  the  tech¬ 
niques  employed  in  this  study  can  be  extended  to  the  problem  of  flow  induced 
vibration,  with  all  non-linearities  geometrical  and  dynamical  included,  of  the 
flexible  aircraft  or  rotating  structures.  From  the  theoretical  point  of  view  since 
this  physical  system  is  a  rich  source  of  static  and  dynamic  instabilities  and 
associated  limit  cycle  motions  it  could  be  used  for  instance  as  a  test  case  for 
assessing  techniques  for  the  study  of  nonlinear  dynamics  and  chaos. 

Subjects  for  the  consideration  in  future  research  include  the  use  of  other  sta¬ 
bility  concepts  such  as  Liapunov  exponents,  and  Poincare  maps.  Application 
of  a  complementary  methods  of  differentiable  dynamics,  in  particular,  of  cen¬ 
ter  manifold  and  bifurcation  theory,  to  analyze  the  problem  &om  a  qualitative 
viewpoint  would  be  helpful.  Knowledge  of  the  generic  structures  of  attract¬ 
ing  sets  in  fV-space,  might  make  the  interpretation  of  numerical  solutions  of 
evolution  equations  considerably  clear. 
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1  Sketch  of  representative  phase  plane  orbits 
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SUMMARY 


The  influence  of  maneuvering  on  the  response  of  a  fluttering  buckled  plate  on  an  aircraft 
has  been  studied.  Assuming  that  the  maneuvering  of  the  aircraft  b  prescribed,  the 
Lagrange  equations  of  motion  for  the  elastic  degree  of  freedom  have  been  derived.  The 
formulation  includes  geometric  non-linearities  associated  with  the  occurrence  of  tensile 
stresses,  as  well  as  coupling  between  the  rigid-body  rotation  of  the  frame  of  reference 
and  the  other  degrees  of  freedom.  This  equation  is  then  tised  to  study  the  response 
of  the  panel  of  an  aircraft  engaged  in  a  pull-up  maneuver.  General  concepts  from  the 
modern  theory  of  dynamical  systems,  with  emphasis  on  chaotic  behavior,  are  presented. 
The  large-amplitude  responses  are  investigated  by  using  the  digital  computer.  Long¬ 
time  histories,  three-dimensional  view  of  orbits,  phase  planes,  and  power  spectra  of 
the  response  are  presented.  We  have  chosen  the  system  parameters  such  that  there 
are  two  equilibrium  positions  for  the  non-maneuvering  case.  The  numerical  simulation 
for  periodic  and  chaotic  response  are  conducted  in  order  to  analyze  the  influence  of 
the  maneuvering  on  the  dynamic  behavior  of  the  panel.  As  the  maneuvering  (load 
factor)  increases,  the  system  exhibits  complicated  dynamic  behavior  including  period- 
multiplying  and  demultiplying  bifurcations  and  chaos. 


SECTION  1 

INTRODUCTORY  REMARKS 


1.1  Overview 

This  report  deals  with  the  response  of  a  flexible  structure,  including  the  dynamical  and 
structural  non-linearities.  The  main  motivation  for  this  work  is  the  aeroelastic  analysis 
of  maneuvering  airplanes;  therefore  aircraft  will  be  emphasized  throughout  this  report. 

The  large  amount  of  material  here  suggested  a  division  of  the  study  into  two  main  areas, 
theoretical  and  numerical.  The  theoretical  approach  is  given  in  Section  2  and  3.  Since 
this  research  is  continuation  of  the  work  by  Morino  and  Baillieul  on  the  geometrically- 
exact  formulation  of  maneuvering-aircraft  dynamics  their  results  are  briefly  discussed 
in  Section  2.  Assuming  that  the  maneuvering  of  the  aircraft  is  prescribed,  and  following 
the  conventional  procedure  of  Lagrangian  mechanics,  the  Lagrange  equations  of  motion 
for  the  elastic  degree  of  freedom  are  derived.  These  are  the  governing  equations  for  the 
physical  system  considered  here.  The  main  ideas  behind  the  modem  theory  of  nonlinear 
systems  are  presented  in  Section  3.  The  reader  tmfamiliar  with  the  subject  hopefully 
will  get  sufficient  information  in  order  for  him  to  be  able  to  follow  the  interpretation 
of  phenomena  appearing  during  the  numerical  experimentation.  The  numerical  results 
are  presented  in  Section  4.  The  governing  equations  are  integrated  numerically  by  using 
a  fourth-order  Runge-Kutta  algorithm.  Some  unusual  non-maneuvering  responses  are 
presented  in  Subsection  4.1.  The  effect  of  varying  the  load  factor  was  studied  in  great 
detail;  only  representative  results  are  presented  in  this  report.  In  order  to  understand 
the  dynamic  problem,  the  character  of  the  solution  has  been  examined  in  Subsection 

4.2  in  physical  terms.  In  Subsection  4.3  the  results  of  the  simulation  of  the  influence  of 
maneuvering  to  a  periodic  response  are  presented.  A  chaotic  responses  are  examined  in 
Subsection  4.4.  Finally,  the  time  dependent  maneuvering  is  considered  in  Subsection 
4.5.  In  examining  and  discussing  the  results,  attention  is  drawn  to  the  routes  to  chaos. 
The  system  exhibits  very  complicated  dynamic  behavior  including  period  multiplying 
and  demultiplying,  bifurcations,  and  chaos.  An  important  result  of  this  work  is  the 
conclusion  that  maneuvering  can  change  the  character  of  the  response  from  periodic, 
or  even  fixed  point,  to  chaotic,  and  vice  verso. 

Although  many  numerical  experiments  were  performed  varying  all  of  the  parameters, 
only  the  results  necessary  to  make  discussion  complete  and  clear  are  presented.  Long 
time  histories,  three-dimensional  view  of  orbits,  phase  planes,  and  power  spectra  of  the 
response  (defined  in  Section  3)  were  the  dynamic  tools  used  in  studying  of  the  system 
considered  here. 
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1.2  Objectives  of  Report 


In  this  report  the  equations  for  the  analysis  of  the  dynamics  of  a  fluttering  buckled 
plate  on  an  aircraft  undergoing  a  pitching  maneuver  have  been  derived.  The  equations 
include  geometric  non-linearities  associated  with  the  occurrence  of  tensile  stresses,  as 
well  as  the  coupling  between  the  rigid-body  rotation  of  the  frame  of  reference  and  the 
the  elastic  degrees  of  freedom.  These  equations  are  then  used  to  study  the  response  of 
the  panel  of  an  aircraft  engaged  in  a  pull-up  maneuver. 

This  research  has  its  own  theoretical  merits  in  understanding  the  influence  of  the  ma¬ 
neuvering  on  the  panel  response  and  occurrence  of  chaos  in  a  system.  The  complexity 
of  the  physical  system  considered  here  justifies  the  use  of  the  simulation  approach. 
Long-time  histories,  three-dimensional  view  of  orbits,  phase  planes,  and  power  spectra 
of  the  response  are  studied.  Finally,  as  a  long  term  goal,  we  believe  that  the  techniques 
employed  in  this  study  can  be  extended  to  the  problem  of  flow-induced  vibration  of 
the  flexible  aircraft,  with  all  the  non-linearities  included  (geometric,  dynamic,  and 
aerodynamic). 

l.S  Review  of  State  of  the  Art 

During  the  last  two  decades,  formulations  for  the  analysis  of  flight  dynamics  and  struc¬ 
tural  dynamics  of  a  flexible  aircraft,  as  well  as  related  formulations  in  the  field  of 
dynamics  of  large  space  structures,  have  been  developed.  An  overview  of  the  recent 
developments  and  trends  in  the  field  of  flexible  aircraft  dynamics  are  given  in  Morino 
and  Baillieul  (1987). 

Since  a  buckled  fluttering  plate  is  a  prototypical  system,  a  brief  review  of  the  recent 
developments  is  in  order.  One  of  the  approaches  to  the  treatment  of  the  flow  induced 
oscillations  relies  on  the  selection  of  relatively  simple  models  which  provides  a  reason¬ 
able  reproduction  of  the  behavior  observed  on  the  real  physical  system.  The  complete 
partial  difierential  equation  is  replaced  by  an  ordinary  differential  equation  on  a  state 
space  of  low  order.  For  example  in  Dowell  and  Pezcshki  (1986,  1987),  Zavodney  and 
Nayfeh  (1988),  or  Zavodney  et.  all  (1989)  Duffing’s  equation  has  been  used  as  a  model 
for  the  sinusoidally  excited  buckled  plate.  In  the  another  approach  the  non-linear  par¬ 
tial  differential  equation  is  recast,  by  Galerkin  or  Rayleigh-Rits  methods,  as  a  set  of 
ordinary  differential  equations  which  are  then  solved  for  specific  initial  conditions  by 
numerical  integration  techniques.  This  approach  is  perhaps  best  illustrated  in  the  work 
of  Dowell  and  Ilgamov  (1988).  A  dissipative  partial  differential  equation  modeling  a 
buckled  beam  was  considered  by  Holmes  and  Marsden  (1961).  They  presented  one  of 
the  very  few  analytical  results  available  on  chaos  in  a  continuous  system. 

Note  that  none  of  the  above  work  has  addressed  the  problem  presented  here:  the 
response  of  a  flexible  plate  of  a  maneuvering  aircraft  including  the  effect  of  the  rigid- 
body  rotation  of  the  frame  of  reference  on  the  deformation. 
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SECTION  2 

THE  EQUATION  OF  MOTION 


The  dynamics  of  a  structure  has  been  formulated  using  Lagrangian  mechanics  in  terms 
of  a  “body  frame  of  reference” .  This  frame  of  reference  is  clearly  defined  in  the  case 
of  a  rigid  body.  In  the  case  of  a  flexible  structure,  the  motion  of  the  structiire  can 
be  expressed  in  terms  of  the  motion  of  the  origin  of  a  body  frame  with  respect  to  a 
chosen  inertial  frame,  the  rigid-body  rotation  of  the  body  frame  around  its  origin,  and 
a  deformation  with  respect  to  a  reference  configuration  rigidly  connected  to  the  body 
frame.  The  orientation  of  the  body  frame  with  respect  to  the  inertial  frame  of  reference 
is  defined  through  the  orthogonal  rotation  matrix. 

We  will  tise  the  following  notations.  An  arrow  above  a  lower  case  letter  designates  a 
vector  in  a  three  dimensional  physical  space.  Superscripted  latin  letters  designate  the 
Cartesian  components  of  vectors  and  tensors.  A  boldface  lower-case  letter  designates  a 
colunm  matrix.  A  boldface  capital  letter  designates  a  matrix;  in  particular,  I  designates 
the  identity  matrix. 

2.1  Lagrange  Equations  of  Motion  for  a  Flexible  Aircraft 

Let  \  s  summarize  the  formulation  and  the  results  presented  in  Morino  and  Baillieul 
(1988);  details  of  the  formulation  are  in  the  original  report.  As  mentioned  before  the 
motion  of  the  aircraft  is  expressed  in  terms  of  the  motion  of  the  origin  of  a  body  frame 
of  reference,  the  rigid-body  rotation  of  the  body  frame  of  reference,  and  a  deformation. 
Indicating  with  x  the  coordinates  in  the  inertial  frame,  with  xq  the  coordinates  of  the 
origin  Pq  of  the  body  frame,  with  s  the  body-axis  coordinates  of  a  point  P  in  the  body 
frame  at  time  t  =  0  (reference  condition)  and  with  u  the  body-axis  components  of  the 
displacement  of  the  structure  with  respect  to  reference  condition,  we  have, 

x(r ,  0  =  xo(0 + 11(0  [»(r ) + u(r .  o]  (2.1) 

where  f**  is  a  system  of  material  coordinates,  (ae.,  a  system  of  coordinates,  in  general 
curvilinear,  that  is  convected  with  the  material  point).  Typically,  these  coincide  nu¬ 
merically  with  the  components  of  a.  Note  that  the  matrix  R  represents  a  rigid-body 
rotation.  Next  let  us  assume  that  a  deformation  is  given  as  a  linear  combination  of 
prescribed  functions,  with  unknown  coefficients,  ae., 

(2-2) 

In  the  rest  of  this  subsection  the  resulting  Lagrangian  equations  of  motion  are  presented 
for  all  the  degrees  of  freedom,  t.e.,  translational  xq,  rotational  R,  and  elastic  «r;  detaUs 
of  the  derivation  are  in  the  original  report. 
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The  following  notation  will  be  used, 

”*  =  (2.3) 

•o  =  (2'^) 

•r  =  J[fyP*rdV  (2.5) 

mr,  =  Jlf^P  4r  =  m„  (2.6) 

b,  =  jjj^pmx^rdV  (2.7) 

b,.  =  J[l^p4r^4.dV  =  -h„  (2.8) 

Jo  =  ffj^p{»^  *!-•%' )dV  =  31  (2.9) 

J.  =  +.^;)]dV  =  J,*  (2.10) 

J.*  =  jjj^p[4!l4,l-\{*r4:  +*Al)\dV  =  3l,=3„  (2.11) 


First,  consider  the  translational  degrees  of  freedom.  Noting  that,  by  definition  of  center 
of  mass,  s  =  m(x(;-xo),  the  Lagrangian  equations  of  motion  for  the  translational  degrees 
of  freedom  are  given,  in  body  axes,  by 


mO(Xo  +  DfS  =  m(vo  +  Ovo)  +  lito  +  n*so 
+  (itir  +  2nu,  +  nu,  +  n*«,)  ■,.  =  Ib  =  ma''  g  +  eo 


(2.12) 


where  fl  =  >t(w)  and  A  denotes  the  operator  that  maps  an  arbitrary  vector  into  a 
skew-symmetric  matrix.  Furthermore,  vo  =  R*  zo,  Dt  denotes  the  operator  !?((•)  = 
Id(-)/dt  +  n(.),  and  eo 

eo=j^td5  (2.13) 

is  the  aerodynamic  resultant. 

Next,  consider  the  rotations'  ^rees  of  fireedom.  With  an  orthogonal-matrix  variable, 
a,  and  following  the  pror  ^  of  Bmllieul  and  Levi  (1087),  the  Lagrange  equation  of 
motion  for  the  rotation  around  the  point  can  be  written  as, 


D|ho  +  s  X  D*vo  =  ms  +  53“'“'  +  ■  x  In 


(2.14) 


This  coincides  with  the  equation  of  conservation  of  the  angular  momentum  obtained 
from  Newtonian  Mechanics.  Note  that,  ho  =  Jw  +  b  is  the  angular  momentum  around 
the  point  Pq  of  the  motion  about  fb,  s&d 


s  X  t  <15 


(2.15) 
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is  the  aerodynamic  moment  (about  the  origin  Pq  of  the  body  frame  of  reference)  acting 
on  the  undeformed  configuration,  whereas 

TOr  =  x  t  dS  (2.16) 

is  the  change  in  moment  due  to  the  change  in  x  caused  by  a  unit  increase  of  the 
Lagrangian  coordinate  u,  (for  a  prescribed  distribution  of  t). 

Finally,  consider  the  Lagrangian  equations  of  motion  for  the  elastic  degrees  of  free¬ 
dom.  Following  the  conventional  procedure  of  Lagrangian  mechanics  one  obtains  the 
equations  of  motion  for  the  Lagrangian  coordinate  tv,  in  body  axes,  as: 

where  £  is  the  clastic  energy,  and  w  are  the  components  of  the  angular  velocity  in  the 
body  frame  of  reference. 

2.2  Dynamics  of  a  Fluttering  Plate  in  a  Maneuvering  Aircraft 

In  this  subsection  the  Lagrange  equations  of  motion  obtained  above  are  specialized  to 
the  case  of  a  fluttering  buckled  plate  on  an  aircraft  tmdergoing  a  pitching  maneuver. 
Since  Wi  want  to  study  the  effect  of  maneuvering  on  the  elastic  response,  we  assume 
the  rigid-body  motion  of  the  frame  of  reference  to  be  prescribed.  Hence  Eqs.  2.12,  2.14 
are  not  used,  and  in  Eq.  2.17  u  is  prescribed  and  therefore  the  inertial  terms  are  linear 
(with  time  dependent  coefficients).  Consider  a  right-handed  body  frame  of  reference 
with  base  vectors  jl,  72173.  We  will  assume  constant  velocity  pull-up  maneuver,  t.e., 
vo  =  -vo7i,  with  u  =  -u;^,  and  vq  =  constant.  Furthermore,  the  modes  4r  tdl  in 
the  direction  73,  and  therefore  the  following  simplifications  are  possible:  nv,  =  MrSrt 
(where  Sr,  designates  the  Kronecker  delta  and  Mr  =  nw  indicates  the  generalized  mass, 
see  Eq.  2.6),  b,,  =  0  (see  Eq.2.8),  and  J,.  =  M,  jr«(I-l3)  (where  I3  =  [jjsjt  see  Eq.  2.11). 
In  addition,  s,  =  »rh*  br  =  Th\is,  the  Lagrangian  equations  for  the  elastic 
degrees  of  freedom  2.17  simplify  considerably  and  are  given  by 

SrWtH)  +  +  5^  nv.fi,  -  53  "’"V.  “•  +  ST  = 

«  « 

with  r  = 

Consider  once  more  a  thin  plate  having  length  a  and  thickness  h  and  undergoing  cylin¬ 
drical  bending  in  response  to  one  side  airflow,  see  Fig.  1.  In  such  a  case  the  axial 
octension,  v,  can  be  written  to  flrst  order  approociznation  as 
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Thus  the  elastic  energy  is 


where  D  is  a  plate  bending  stiffness,  and  K  a  spring  constant  per  unit  spanwise  length 
of  panel.  Substituting  Eq.  2.2  into  the  elastic  energy,  Eq.  2.20,  and  differentiating 
with  respect  to  u,  one  obtains 

^  =  (2-21) 
«  m  • 

where 

=  rCi^rdz  =  r  4^y;dz.  (2.22) 

Jo  Jo 

In  the  last  equation,  an  integration  by  parts  has  been  used  with  =  0  at  z  =  0. 

We  will  assume  that  the  plate  is  exposed  to  a  in-plane  tensile  load,  JV,,  to  a  static  pres- 
stire  difference  across  the  plate,  AP,  and  that  the  exciting  dynamic  pressure  difference 
is  given  by  the  quasi-steady  aerodynamic  theory,  t.e., 


P-Poo 


2q  fdtt  Af*-2  1  du 
0  dz'^  -I  vq  dt 


(2.23) 


where  q  =  pvll 2  the  dynamic  pressure,  M  is  the  Mach  number,  and  =  VAf*  -  1, 
Bisplinghoff  and  Ashley  (1962). 

Following  the  conventional  procedure  of  Lagrangian  mechanics,  t.e.,  calculating  the 
virtual  work  done  by  the  external  forces,  combining  these  expressions  with  Eqs.  2.2, 
2.18,  and  2.21  one  obtains 


where 


»rV3Vo  +  fc,tD  +  5^  ^  H 

•  #  9  m  9 


,  P 


53(««Pr.  +  OrA-.) 


APjv, 


(2.24) 


Prs  = 

l\s4>rdx, 

JO 

PT  = 

Jo 

Pr,  = 

L 

fVj  = 

-  r  4>‘l4rdz  =  f 

(2.25) 

In  the  last  equation  an  integration  by  parts  has  been  used  with  ^  =  0  at  z  =  0. 

In  order  to  put  Eq.  2.24  in  a  dimensionless  form,  we  introduce  the  following  dimen¬ 
sionless  parameters  and  coordinates: 

_«r  D 

h'  [pmha*)  ’ 

O  -  1  - 

”  ~  tTo  ’  /3D  ’ 

iV.o* 


M  = 


pa 

Pn.h' 


a  = 


Dh 

Ka 


{Ka  +  Ek)' 


(2.26) 


Substituting  4>r  =  •ixi{ritz/a)  into  Eq.  2.25,  and  the  result  in  Eq.  2.24,  the  equation  of 
motion  in  dimensionless  matrix  form  rends 


ik  +  Sii+  Ga  +  nd  +  Ov  =  p  -I-  f(o)  +  7n*a 
where  G  =  [y,„],  d  =  {d,},  v  =  {»,},  p  =  {p,},  f  =  {/,},  and 
fr  =  3(rx)<a(l  - 


(2.27) 


Srm  =  2A 


rm 


r2- 


(1  -  (-!)'♦”)  +  ((nr)»Ji.  +  (nr)‘X™. 


0x .  [1  -  (-m 

p,=2P^— 


(2.28) 


with  m,  r  =  1, ...,  N  for  the  elastic  degrees  of  freedom.  Eqs.  2.27  are  a  set  of  ordinary, 
nonlinear  differential  equations  in  time.  The  cubic  type  of  non-linearities  are  of  a 
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geometric  origin  associated  with  the  occurrence  of  tensile  stresses  in  the  middle  surface. 
The  coupling  between  the  rigid-body  rotation  of  the  frame  of  reference  and  the  elastic 
degrees  of  freedom  is  represented  by  the  fourth  term.  Note  that  by  mtmiTniTig  0  =  0 
Eqs.  2.27  properly  reduce  to  a  classical  panel  flutter  equations,  see  Bolotin  1963,  or 
Dowell  1966,  for  example.  This  set  of  equations  have  been  solved  by  a  direct  numerical 
integration. 
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SECTION  3 

ANALYSIS  OF  NONLINEAR  SYSTEMS 

In  this  section  we  introduce  the  main  ideas  behind  the  modem  theory  of  nonlinear 
systems  or,  more  precisely,  the  qualitative  theory  of  differentiable  dynamical  systems. 
Our  goal  is  to  provide  sufficient  information  to  the  reader  who  is  not  familiar  with  a 
subject  in  order  for  him  to  be  able  to  follow  the  interpretation  of  phenomena  appearing 
during  the  numerical  experimentation.  Readers  desiring  more  detailed  discussion  in 
nonlinear  dynamical  systems  should  consult  books  such  as  Berge,  Pomeau,  and  Vidal 
(1984),  Holden  (1986),  Moon  (1987).  Readers  looking  for  the  more  rigoroiis  formula¬ 
tions  should  consult  books  such  Guckenheimer  and  Holmes  (1983),  Anosov  and  Arnold 
(1988),  or  Arnold  (1989). 

It  should  be  emphasized  that,  during  the  last  decade  there  has  been  an  increased 
interest  in  the  mathematics  of  chaotic  dynamics.  The  literature  is  quite  extensive  and 
a  review  is  not  attempted  here.  However,  the  book  on  chaos  by  Hao  Bai-Lin  (1984) 
should  be  mentioned:  it  includes  an  introduction  as  well  as  reprints  of  several  important 
publications  on  chaos  (forty one  papers),  and  over  five  hundred  references. 

5.1  Classical  Nonlinear  Dynamics  Theory 

In  this  subsection,  a  brief  review  of  nonlinear  vibration  theory  is  presented.  The  idea 
is  simply  to  define  and  review  the  basic  concepts  of  classical  nonlinear  vibration,  so 
that  we  may  later  be  able  to  contrast  these  with  the  concepts  of  the  modem  theory  of 
nonlinear  systems. 

5.1.1  Vector  Fields,  Differential  Equations 

Consider  a  real,  finite  dimensional,  linear  space  V  =  JT*.  A  vector  field,  defined  in  a 
domain  U  of  the  space  V,  is  a  map  which  associates  to  each  point  x  £  U  h.  vector  v, 
based  at  x,  in  the  space  V.  The  equation 

x  =  v(x,0,  X£UCV,  (3.1) 

is  called  the  differential  equation  corresponding  to  the  vector  field  v.  In  Eq.  3.2  the 
following  notation  has  been  used  (')  =  d(  )/cft,  and  the  variable  (  is  called  '^ime”.  If 
the  right  hand  side  is  time  independoit 

i  =  ▼(*),  (3.2) 

the  differential  equation  is  called  autonomous^. 

'Note  that  tha  tanas  ^toooaiogs”  aad  "Ubm*  fai  tha  abeva  daflaithaw  cams  «at  of  phftics;  for  iaataaca,  tha 
aaolotioa  law  tit  a  idijaically  aatoBomona  srstam  ((.a.,  a  ajatam  which  doas  aot  lataract  with  othar  srataais), 
doas  aot  aanall  j  dspaad  oa  thna 
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The  domain  U  is  called  the  phase  space  of  this  equation,  and  x(t)  is  the  state  vector. 
A  system  of  differential  equations  such  as  3.1  or  3.2  is  called  a  flow  in  the  domain  U. 

Often  we  are  given  an  initial  condition  x(0)  =  xo£U  and  we  seek  a  solution  x(xo,  t)  or 
simply  x(t)  which  satisfies  3.1  or  3.2.  An  integral  curve  of  a  differential  equation  is  the 
graph  of  a  solution;  a  phase  curve  or  tr^ectory  or  orbit  is  the  projection  of  an  integral 
curve  on  the  phase  space  along  the  t-axis. 

One  wishes  to  analyze  the  evolution  of  arbitrary  physical  systems  described  by  the 
differential  equations  3.1,  or  3.2.  The  objective  of  the  analysis  is  to  characterize  in 
broad  outline  the  critical  aspects  of  system  behavior,  not  the  details.  These  critical 
aspects  are  discussed  below. 

S.1.2  Dissipative  Systems 

Let  V  be  a  vector  field.  If  div  v  =  0,  then  the  phase  flow  of  Eq.3.2  preserves  volume.  If 
div  V  <  0  everywhere  in  the  phase  space,  then  the  phase  space  flow  decreases  volume 
as  time  increases  (t.e.,  the  volume  corresponding  to  a  set  of  points  at  t  =  0  decreases  it 
time).  L.  this  case  the  equation  is  said  to  be  dissipative. 

Note  that  contraction  of  the  volume  can  be  obtained  by  reduction  of  all  lengths,  or 
by  a  decrease  of  some  of  the  lengths,  accompanied  by  the  less  rapid  increase  of  other 
lengths.  This  remark  is  important  for  showing  that  divergence  of  phase  trid^ctories  in 
some  direction  remains  possible,  even  in  a  dissipative  system.  Knowledge  of  the  way 
in  which  the  decreasing  of  the  volume  takes  place,  and  of  the  speed  at  which  it  occurs 
are  essential  for  a  complete  dynamical  description. 

S.l.S  Singular  Points 

The  concept  of  an  singular  point,  which  was  used  extensively  in  linear  dynamic  systems, 
carries  over  directly  to  nonlinear  dynamic  systems. 

A  singular  point  of  a  vector  field  v  is  a  point  at  which  the  vector  field  vanishes.  A  sin¬ 
gular  point  of  a  differential  equations  3.1,  or  3.2  is  a  singular  point  of  the  corresponding 
vector  field  v. 

The  singular  points  of  a  differential  equation  are  sometimes  called  equilibrium  points. 

Since  equilibrium  points  are  solutions  to  nonlinear  equations,  finding  such  solutions 
is  somewhat  more  of  an  accomplishment  than  for  the  linear  case.  Furthemoore,  the 
equilibrium  points  distribution  may  be  more  complex  than  for  a  linear  case;  the  system 
may  have  none,  one,  or  more  equilibrium  points  and  they  may  be  in  any  pattern  in 
phase  space. 

Ultimately,  however,  as  in  the  linear  case,  interest  centers  not  just  on  the  existence  of 
equilibria,  but  also  on  their  stability  properties. 

The  problem  of  stability  depends  on  whether  the  system  is  dissipative  or  conservative. 
We  will  deal  with  the  dissipative  systems  only. 
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5.1.4  Stabmty 

Let  xo  be  a  singular  point  of  a  vector  field  v.  A  stationary  solution  of  an  autonomous 
differential  equation  3.2  (t.e.,  a  solution  whose  image  is  an  equilibrium  point)  is  said 
to  be  Lyapunov  stable  if  all  solutions  of  the  equation,  with  initial  conditions  in  a  suf¬ 
ficiently  small  neighborhood  of  the  equilibrium  point,  are  defined  for  all  t  >  0  and 
converge  uniformly  (with  respect  to  time)  to  the  stationary  solution  as  the  initial  con¬ 
ditions  tend  to  the  equilibrium  point.  Fig.  2.a. 

A  stationary  solution  is  said  to  be  asymptotically  stable  if  it  is  Lyapimov  stable  and 
if,  in  addition,  all  solutions  with  initial  conditions  sufficiently  close  to  the  equilibrium 
point  under  consideration  tend  to  this  equilibrium  point  as  f  -*  oo. 

A  stationary  solution  is  marginally  stable  if  it  is  stable  but  not  asymptotically  stable. 

A  solution  is  unstable  if  it  is  not  stable. 

Lyapunov  stability  of  any  solution  of  any  differential  equation  (autonomous  or  not)  is 
defined  in  similar  manner.  That  is,  the  solutions  for  all  t  >  0  are  required  to  converge 
uniformly  to  the  solution  in  question  as  their  initial  values  at  (  =  0  tend  to  the  initial 
value  of  the  solution  in  question. 

According  to  the  basic  definitions,  stability  of  a  stationary  solution  is  a  local  property 
of  the  vector  field  (defining  the  differential  equation)  at  the  corresponding  equilibrium 
point.  Therefore,  to  conduct  an  analyse  of  stability  it  is  often  theoretically  legitimate, 
and  mathematically  convenient  to  replace  full  nonlinear  description  3.2  with  the  linear 
approximation  at  the  singular  point. 

5.1.5  Linearisation  and  Stability 

Let  Xo  be  a  singular  point  of  a  differentiable  vector  field  v,  and  let  A  be  the  Jacobian 
matrix  of  the  vector  field  v 


The  equation 


A  —  [oij],  Oij  — 


X  =  Ax, 


is  called  the  linearization  of  Eq.  3.2  at  the  singular  point  xo.  A  singular  point  of  a 
vector  field  is  said  to  be  non-degenerate  if  the  Jacobian  matrix  is  nonsingular  at  this 
point.  The  field  Ax  is  called  the  linear  part  of  the  field  t  at  xo. 

If  all  the  eigenvalues  of  the  linear  part  of  a  vector  field  v  at  a  singular  point  have 
negative  real  part,  then  the  singular  point  xo  is  asymptotically  stable  for  the  nonlinear 
system,  ff  at  least  one  eigenvalue  has  a  positive  real  part,  the  singular  point  is  not 
Lyapunov  stable,  ff  any  one  of  the  eigenvalues  has  sero  real  part,  then  stability  cannot 
be  determined  by  linearization. 


We  say  that  a  singular  point  is  hyperbolic  if  the  linear  part  of  the  vector  field  v  has 
no  eigenvalue  with  real  part  zero.  A  hyperbolic  singular  point  is  either  unstable  or 
asymptotically  stable. 

If  all  eigenvalues  of  the  linear  part  of  a  vector  field  v  at  a  hyperbolic  singular  point 
have  negative  real  parts,  the  singular  point  is  called  a  sink;  if  they  are  all  positive  the 
singular  point  is  a  source.  If  both  signs  occur,  the  singular  point  is  a  saddle  point. 

S.1.6  Cycles 

It  is  well  known  that  in  a  physical  system  whose  law  of  evolution  does  not  change 
in  time,  periodic  regimes  may  be  established.  A  mathematical  description  of  these 
phenomenon  is  given  by  the  theory  of  cycles  developed  by  H.  Poincari. 

A  closed  phase  curve  of  some  fiow  in  If,  arising  from  a  nonlinear  vector  field  v,  is 
called  a  cycle.  Choose  a  point  A  on  the  ctirve  and  draw  a  local  cross  section  Eel/, 
of  dimension  n  -  1,  where  n  denotes  the  dimension  of  the  phase  space.  Fig.  2.b.  The 
hypersurface,  £  must  be  chosen  so  that  the  flow  is  everywhere  transverse  to  it,  t.e., 
v(x)  •  n(x)  ^  0  for  all  x  €  £,  where  n(x)  is  the  unit  normal  to  £  at  x.  Phase  curves 
which  start  at  points  of  £  sufficiently  close  to  the  original  point  of  the  cycle  return  to 
£.  Thus  we  obtain  a  set  of  points  around  A  corresponding  to  the  starting  points.  This 
relationship  defines  a  map  of  £  into  itself  which  is  called  first  return  map  or  Poincare 
map. 

Th^  mathematical  study  of  such  maps  is  similar  to  that  for  difTerentia)  equations.  One 
can  find  equilibrium  points  of  the  map  and  one  can  classify  these  equilibrium  points 
by  the  study  of  linearized  maps  about  the  equilibrium  points. 

It  is  apparent  that  the  point  A  on  the  cycle  is  an  equilibrium  point  of  the  Poincari 
map,  and  that  the  stability  of  an  equilibrium  point  for  the  Poincari  map  reflects  the 
stability  of  closed  phase  curve  for  the  flow.  Furthermore,  Poincari  maps  have  the  same 
kind  of  topological  properties  as  the  flow  from  which  they  arise.  For  example  if  the 
flow,  Eq.  3.2,  is  dissipative,  so  that  volumes  in  the  phase  space  are  contracted,  then 
the  map  contracts  areas  in  the  plane  £. 

A  closed  phase  ctirve  is  called  a  limit  cycle  if  the  flxed  point  of  the  corresponding  first 
return  map  is  isolated. 

Multipliers  of  a  cycle  are  defined  to  be  the  eigenvalues  of  the  linear  part  of  the  first 
return  map  at  the  fixed  point  corresponding  to  the  cycle. 

A  cycle  is  said  to  be  Lyapunov  (orbitally)  if,  for  an  arbitrarily  small  neighborhood  C  of 
the  cycle,  all  tr^ectories  which  start  in  a  sufliciently  small  neighborhood  of  the  cycle 
do  not  emerge  from  C. 

A  cycle  is  said  to  be  Lyapunov  asymptotically  stable  if  it  is  Lyapunov  stable,  and  if  all 
the  phase  curves  with  initial  condition  sufficiently  close  to  the  C3rele  approach  the  cycle 
asymptotically  as  (  -»  oo. 
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If  the  modulxis  of  every  multiplier  of  a  cycle  is  less  then  one,  then  the  cycle  is  Lyapunov 
asymptotically  stable. 

S.2  Chaos  in  Dissipative  Systems 

In  this  subsection  we  will  define  the  different  types  of  time-evolution  of  a  dissipative 
systems,  including  those  which  at  first  seem  totally  disordered,  chaotic. 

5.2.1  Asymptotic  Behavior  of  a  Dissipative  System,  Attractors 

Consider  an  autonomous  dissipative  system  in  the  phase  space  U.  By  definition  of 
dissipative  system,  the  phase  space  fiow  decreases  volume  as  time  increases.  Continued 
decreasing  of  volume  in  phase  space  means  that  as  f  -»  oo  the  phase  curves  of  an  evo¬ 
lutionary  process  get  closer  and  closer  to  a  compact  aet  M  cU,  the  attractor.  Motions 
starting  in  some  volume  of  phase  space,  the  basin  of  attraction  for  the  attractor,  can 
follow  complicated  transients  but,  as  t  -»  oo,  they  finally  approach  the  attractor.  The 
attractor  M  is  invariant  tmder  the  fiow;  any  motion  on  the  attractor  is  confined  to  the 
attractor.  There  is  not  transient  motion  on  the  attractor;  all  transients  in  a  dissipative 
system  are  in  the  basin  of  attraction  of  an  attractor  on  the  approach  to  an  attractor, 
and  are  not  on  attractors. 

Asymptotically  stable  singular  points  and  Lyapunov  asymptotically  stable  cycles  are 
particular  cases  of  attractors. 

5.2. 2  Strange  Atvractors 

A  strange  attractor  is  an  attracting  set  that  is  not  a  singular  point  nor  a  periodic  cycle, 
and  for  which  almost  all  solutions  display  sensitive  dependence  on  initial  conditions. 
Thus  a  triy^^o^  asymptotic  to  a  strange  attractor  M  rapidly  approaches  the  subset 
M  cU  but  thereafter  wanders  erratically  about  the  "surface’*  of  M. 

S.2.S  Characterising  the  Attractors 

There  are  several  ways  to  identify  and  to  characterize  the  nature  of  a  dynamical  regime. 
We  will  expl^  here  the  power  spectrum  method  used  latter  during  numerical  experi¬ 
mentation. 

Let  ▼  be  a  nonlinear  vector  field  in  the  domain  U.  Let  s(t)  be  a  component  of  the  state 
vector  X.  Note  that  numerical  simulation  generally  provides  a  discrete  sequence  of  real 
nuxnbers  j  €  Z  regularly  spaced  at  time  intervals  At.  In  practice  this  sequence  of 
numbers  is  necessarily  finite;  suppose  that  we  have  N  consecutive  sampled  values.  We 
define  the  Fourier  transform  of  a  discrete  time  series  {cj}  to  be  the  operation  creating 
a  corresponding  discrete  series  {X{wi,)}  such  that 

=  (3-5) 
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over  discrete  set  of  frequencies 


Wfc  =  k  =  -N/2, ....  N/2.  (3.6) 

The  power  spectrum  of  Xj  for  j  =  is  defined  as 

PuM  =  K\XM\\  (3.7) 

where  K  is  a  scale  factor,  see  Priestley  (1981)  for  details,^. 

Having  defined  the  discrete  Fourier  transform  and  power  spectrum  of  the  discrete  series, 
we  return  to  the  explanation  of  a  relationship  between  the  nature  of  the  series  and  the 
form  of  its  power  spectrum. 

In  essence,  the  spectrum  of  a  periodic  series  of  period  T  is  made  up  of  a  peak  at  the 
frequency  1/T.  A  quasiperiodic  sequence  has  r  fundamental  frequencies  Wj,  j  =  1,  ...,r; 
its  Fourier  spectrum  contains  components  at  all  frequencies  of  the  form 


jwimi  +  (021x^2  +  ...  +  w,m,|,  (3.8) 

where  my  are  arbitrary  integers.  When  the  series  is  neither  periodic  nor  quasiperiodic, 
it  is  called  aperiodic,  the  Fourier  spectrum  is  then  continuous. 

S.2.4  Pseudo-Definition  of  Chaos 

Given  that  no  precise  scientific  defin  .tion  exists  for  the  noun  "chaos”  or  for  the  adjective 
"chaotic”,  we  will  consider  these  words  to  be  synonymous  with  certain  properties. 
Following  Berge,  Pomeau,  and  Vidal  (1984),  we  will  say  that  the  dynamical  regime  is 
chaotic  if  its  power  spectrum  contains  a  continuous  part-a  broad  band-regardless  of 
the  possible  presence  of  the  peaks. 

The  real  difficulty  is  that  a  power  spectrum  that  looks  continuous  cannot  be  auto¬ 
matically  attributed  to  an  chaotic  regime  because  this  is  also  the  appearance  of  the 
spectrum  of  a  quasiperiodic  regime  with  a  very  high  number  of  frequencies. 

The  definition  is  pragmatic;  it  lacks  rigor  and  contains  unavoidable  ambiguities. 

S.S  Bifurcation  Theory- Ways  in  Which  Chaos  Appears 

So  far  we  have  been  dealing  with  an  individual  systems,  and  now  let  us  consider  systems 
whose  behavior,  expressed  in  terms  state  vector  x  E  is  in  some  way  controlled  by 
external  variables.  The  latter  are  collected  in  a  control  parameter  p;  one  can  often  take 
p  €  JP”.  The  family  of  differential  equations 

X  =  vp,  X  6  C  JT,  p  €  ir,  (3.9) 

*Net«  that  Um  fut  FonrMr  transfonn  alcoritlim  haa  baaa  oaad  to  compata  tha  diacrata  Fooilar  tranafonn  doring 
tha  Boinarical  aimnlatioii. 
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generates  a  flow  or  vector  field  on  U. 

The  state  space  in  this  case  can  be  pictured  as  filled  with  flow  lines  defined  by  Eq.  3.9 
which  approaches  the  attractors  in  U  as  t  -»  oo;  there  may  also  be  saddle  type  solutions 
in  V.  The  structure  of  the  collection  of  attracting  sets  A  in  U  is  clearly  important; 
one  may  wish  to  know  what  qualitative  solution  types  can  exists  and  in  particular  to 
classify  the  ways  in  which  A  might  change  with  /x  (as  v^  changes).  Qualitative  changes 
in  Yfi  are  referred  to  as  bifurcations. 

S.3.1  Structural  Stability 

It  may  be  helpful,  before  discrissing  bifurcations,  to  give  a  working  definition  of  struc¬ 
tural  stability  and  codimension. 

A  vector  field  is  structurally  stable  if  its  qualitative  structure  is  not  destroyed  by  a 
small  perturbation  in  v^. 

It  is  clear  that  a  vector  field  possessing  a  non-hyperbolic  fixed  point  cannot  be  struc¬ 
turally  stable,  since  a  small  perturbation  can  remove  it  (if  the  linearized  matrix  is 
noninvertible,  having  a  zero  eigenvalue),  or  turn  it  into  a  hyperbolic  sink,  a  saddle,  or 
a  source  (if  the  matrix  has  purely  imaginary  eigenvalues).  Similar  observations  apply 
to  periodic  orbits  and  we  conclude:  if  the  flow  is  structurally  stable  all  fixed  points 
and  closed  orbits  must  be  hyperbolic.  This  is  a  necessary  but  not  sufficient  condition 
to  guarantee  structural  stability,  since  some  global  effects  also  come  into  play.  For 
example,  in  the  case  of  a  structurally  stable  systeri  on  the  plane,  besides  the  above 
conditions  there  must  also  be  no  orbits  connecting  saddle  points. 

The  codimension  of  a  bifurcation  will  be  the  smallest  dimension  of  a  parameter  space 
which  contains  the  bifurcation  in  a  persistent  way. 

Let  us  illustrate  this  concept  in  a  three-dimensional  parameter  space  (m  =  3).  The 
lociis  of  points  in  parameter  space  where  some  condition  is  satisfied,  say  change  of 
stability  of  a  solution,  is  generally  a  surface  S.  Almost  any  line  D  (dimension  one) 
crosses  this  surface;  by  moving  along  17,  we  eventually  intersect  S  and  observe  the 
corresponding  bifurcation.  Thus  we  say  that  the  bifurcation  is  of  codimension  one.  If 
we  have  two  conditions,  these  will  both  be  verified  only  along  the  line  of  intersection  L 
of  the  two  surfaces  £i  and  £3  associated  with  the  two  conditions.  An  arbitrary  line  has 
no  chance  of  intersecting  X;  we  must  move  along  a  surface  (dimension  two)  to  cross  a 
line  of  bifurcations.  This  is  an  example  of  the  codimension  two  bifurcation. 

What  we  have  just  said  for  m  =  3  can  be  generalized  to  any  other  value  of  m  simply  by 
replacing  the  words:  space,  surface,  etc.,  by  hyperspace,  hypersurface,  etc.  In  a  sense 
the  codimension  of  a  degenerate  fixed  point  (limit  cycle,  etc.)  describes  its  degree 
of  structural  stability;  hyperbolic  fixed  points  are  of  codimension  0.  The  higher  the 
codimension  the  more  varied  and  complex  are  the  ways  in  which  the  singularity  can 
unfold  in  response  to  perturbations.  To  generically  capture  a  codimension  it  singularity 
one  needs  a  k-parameter  vector  field. 
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S.S.2  Bifurcation 

A  Bifurcation  occurs  when  a  vector  field  passes  through  a  regime  of  structural  instability 
and  suffers  a  qualitative  change.  The  parameter  values  at  which  this  change  occur  are 
called  bifurcation  values. 

There  are  two  main  approaches  open  to  us,  local  and  global  analysis.  Local  analysis  is 
generally  performed  by  studying  the  vector  field  near  the  bifurcating  equilibrium  point 
or  closed  orbit,  and  bifurcating  solutions  are  also  found  in  a  neighborhood  of  that  limit 
set.  If  the  dynamical  properties  involve  global  aspects  of  fiows,  they  cannot  be  deduced 
from  local  information,  and  a  global  analysis  is  necessary.  Note  that  even  a  study  of 
local  two-parameter  bifurcations  requires  an  understanding  of  global  bifurcations,  since 
they  occur  naturally  in  two-parameter  families.  Also  note  that  the  global  combinations 
of  local  events  may  be  extremely  complex. 

We  will  examine  only  local  bifurcations  of  individual  equilibria  and  periodic  orbits 
referring  the  reader  interested  in  bifurcations  of  codimension  two  and  codimension 
three  to  more  specialized  articles  and  books. 

Codiznension  One  Bifurcations  of  Equilibria 
Consider  a  fiow 

X  =  v^(x),  (3.10) 

on  the  smooth  vector  field  v^,  with  x  €  €  R.  Asstime  that  at  =  /lo,  x  =  xo,  Eq. 

3.10  has  an  equilibrium  point  at  which  there  is  a  zero  eigenvalue  (for  the  linear  part  of 
the  vector  field).  It  can  be  shown  (Tie  Center  Manifold  Theorem’  )  that  the  dynamic 
on  the  centre  manifold  is  locally  equivalent  to  one  in  which  c  is  one  dimensional.  The 
possible  bifurcations  are  represented  by  the  following  three  differential  equations  which 
depend  on  a  single  parameter  fi  (see  Table  1.). 

Table  1.  The  local  models  on  the  centre  manifolds 


The  local  model 

Name  of  bifarcatioa 

i  =  M  d:  ** 

Saddle-aode  (or  tanuas  poiat) 

Ikaaacritical 

i  —  It9±  a* 

Pitchfork 

Taking  the  minus  sign  in  the  case  of  saddle-node  bifurcation  one  immediately  sees 
that  for  /t  <  0  there  are  no  fixed  points  and  that  for  >  0  there  are  two,  a  saddle  at 
(-  y^,0)and  a  sink  at  (+  ^,0).  Similar  analysis  is  possible  for  other  cases,  and  the 
evolution  of  the  attractors  may  be  pictured  as  in  Figs.  3.a,  3.b,  and  3.c. 

Consider  now  a  system  given  by  Eq.  3.10  with  a  parameter  value  ito  end  equilibrium 
xq  at  which  the  linear  part  of  vector  field  v,,  has  a  simple  pair  of  pure  imaginary 

*W«  hsT*  aot  dafinad  or  dianiaaad  tha  cantar  thaaiam  bacamaa  It  ia  aat  rdavaat  for  tha  diaauaioiu 

latar  ia  thk  laport  othar  than  to  oadarataad  tha  tjpical  bitecatioaa  aacoBatand.  lataiaatad  raadan  ahonld 
coaaalt  Gackaabaiaiar  aad  Holmaa  (19S3),  for  iatrodoctiooa  to  tha  aabjaet. 
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eigenvalues  and  no  other  eigenvalues  with  zeto  real  parts.  Here  the  local  model  on  the 
centre  manifold  is  given  by 

ii  =  Z2  ±  zi[n  -  {xl  +  xl)],  (3.11) 

ij  =  -*i  ±*2[/i  -  (*J  +  *1)]. 

Taking  the  pliis  sign  one  can  verify  that  there  is  a  spiral  sink  at  (0,0)  for  /i  <  0  and  a 
soxirce  at  (0,0)  surrounded  by  a  limit  cycle  for  n  >  0.  The  limit  cycle  evolves  contin¬ 
uously  from  the  centre  at  (0,0)  for  /t  =  0  as  increases.  Fig.  3.d.  The  bifurcation  is 
named  after  E.  Hopf,  who  first  classified  it. 

Note  that,  in  all  the  cases  considered,  the  nonlinear  terms  have  had  an  effect  opposite 
to  that  of  the  instability  caused  by  the  lower  order  terms.  The  bifurcation  in  this 
case  is  called  supercritical  (also  normal).  However,  if  the  lowest-order  nonlinear  terms 
also  have  a  destabilizing  infiuence  on  the  solution  (alternative  signs  in  the  expressions 
listed  above),  the  bifucation  is  subcritical  or  inverse.  We  then  obtain  the  bifurcation 
diagrams  presented  in  Fig.  4. 

Codimension  One  Bifurcations  of  Periodic  Orbits 

The  usual  approach  consists  in  computing  Poincari  return  maps  and  then  in  applying 
the  results  of  the  previous  subsection  for  these  discrete  dynamical  systems. 

There  are  three  ways  in  which  a  fixed  point  of  Poincari  return  maps  may  fail  to  be 
hyperbolic:  at  least  one  eigenvalue  of  the  linear  part  of  the  map  at  the  fixed  point  is 
+1,  -1,  or  a  pair  of  complex  eigenvalues  A,  A,  with  |A1  =  1. 

The  bifurcation  theory  for  fixed  points  with  eigenvalue  1  does  not  change  very  much 
from  what  has  been  explained  for  equilibria  with  eigenvalues  0.  The  one  parameter  fam¬ 
ily  has  a  two-dimensional  center  manifold  on  which  the  dynamics  is  locally  equivalent 
to  either  the  saddle-node,  transcritical  or  pitchfork  bifurcation. 

Bifurcations  with  eigenvalue  -1  do  not  have  an  analogue  for  equilibria,  and  are  as¬ 
sociated  with  flip  bifurcations,  also  referred  to  as  period  doubling  or  subharmonic 
bifurcation.  Let  us  see  what  happens  when  an  eigenvalue  approaches  the  critical  value 
-1  and  Anally  attains  it.  Let  7  be  a  limit  cycle  of  period  T,  and  let  O  be  a  flxed  point 
of  the  corresponding  flrst  return  map.  Before  the  bifurcation,  any  displacement  from 
O  say  zo  decreases  with  every  period,  since  its  value  x  is  multiplied  by  the  negative 
factor  of  absolute  value  smaller  than  one.  The  limit  cycle  is  therefore  linearly  stable. 
When  the  eigenvalues  becomes  equal  to  -1,  the  situations  change.  The  modulus  of 
any  initial  displacement  xq  is  conserved  while  its  sign  changes  by  each  intersection.  We 
therefore  see  a  new  periodic  orbit,  whose  period  is  twice  that  of  the  original  orbit.  Note 
that  conditions  for  the  stability  of  this  period  two  orbits  will  not  be  considered  here; 
interested  readers  should  consult  Guckenheimer  and  Holmes  (1083)  for  detauls. 

We  now  turn  to  the  bifurcations  of  a  periodic  orbit  at  which  there  are  complex  eigen¬ 
values  A,X  with  |A|  =  1.  Analogy  wi^  the  theory  of  the  Hopf  bifurcation  suggested 
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the  name  Hopf  bifurcation,  despite  several  features  which  distinguish  the  two  bifurca¬ 
tions.  The  fixed  point  of  the  first  return  map  is  replaced  by  a  set  of  points  on  a  circle. 
The  fiows  near  the  bifurcation  have  quasiperiodic  orbits,  and  more  subtle  analysis  is 
required  to  capture  these,  see  Guckenheimer  and  Holmes  (1983)  for  details. 

Globalisation 

Note  that  the  essential  models  of  biftircation  outlined  above  are  only  valid  locally,  t.e., 
in  some  neighborhood  F  C  U  of  the  bifurcating  fixed  point  and  in  some  neighborhood 
P  c  of  the  bifurcational  value  of  the  control  parameter  /i.  If  a  number  of  distinct 
bifurcations  occur,  then  one  must  consider  problem  of  fitting  the  various  local  modeb 
together  into  an  overall  scheme.  Thb  problem  of  globalization  b  largely  open  and 
appears  generally  diflicult.  For  more  information  readers  are  referred  to  Arnold  (1989) 
or  Guckenheimer  and  Holmes  (1983). 

S.S.S  Routes  to  Chaos 

One  of  the  main  interests  of  current  theoretical  research  in  nonlinear  dynamics  b  in 
characterizing  the  nature  of  the  transition  from  regular  periodic  dynamics  to  irregular, 
chaotic  dynamics.  The  structure  of  the  collection  of  attracting  sets  and  in  particular 
the  ways  in  which  they  might  change  with  changing  the  parameter  n  b  difficult  and 
largely  open  problem.  In  the  case  of  one-parameter  systems  we  may  quote  Arnold 
(1986)  pp.  24, 

“The  transition  from  a  stable  equilibrium  state  to  a  strange  attractor  can 
take  place  both  by  means  of  a  jump,  or  after  a  mild  loss  of  stability,  Fig.  5. 

In  the  latter  case  the  stable  cycle  which  was  created  itself  loses  its  stability. 

The  loss  of  stability  of  a  cycle  in  a  generic  one-parameter  family  of  systems 
can  take  place  in  a  number  of  ways:  1)  collbion  with  an  unstable  cycle,  Fig.6, 

2)  doubling,  Fig.7,  and  3)  “the  birth  or  death  of  a  torus,”  Fig.  8. 

The  proof  of  the  exbtence  of  such  a  routes  does  not  follow  from  a  rigorously  demon¬ 
strated  theory,  but  rather  from  a  combination  of  guesses  and  numerical  experimenta¬ 
tion. 
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SECTION  4 

NUMERICAL  EXPERIMENTATION 


In  this  report  nvimerical  simulation  for  periodic  and  chaotic  response  are  conducted  in 
order  to  analyze  a  behavior  of  a  fluttering  panel  of  a  maneuvering  airplane.  Equations 
2.27  were  simulated  on  the  digital  computer  by  using  a  fourth  order  Rimge-Kutta 
algorithm.  For  all  of  results  reported  below,  njM  =  0.023, 1/  =  0.3,  and  the  plate 
response  is  calculated  at  z/a  =  0.75. 

4.1  Response  Without  Maneuvering 

Extensive  numerical  experiments  were  performed  by  varying  all  the  parameters,  for 
maneuvering  as  well  as  non-maneuvering  conditions.  Although  our  main  interest  is  to 
study  the  influence  of  maneuvering  on  the  panel  response,  here  we  will  present  some 
unusual  non-maneuvering  responses. 

According  to  Dowell  (1088)  the  distinctive  types  of  motion  which  may  occur  in  a 
phase  plane  (u,  u)  (with  u  being  the  panel  vertical  displacement  at  a  given  point,  e.^., 
z/c  =  0.75)  are  the  following,  see  Fig.  9:  (l)  For  a  small  in-plane  force,  A.,  a  flow 
velocity  might  not  be  sufflciently  large  to  disturb  the  zero  equilibrium  position.  Fig. 
O.a,  (2)  The  buckled  plate  corresponds  to  two  points  in  the  phase  plane.  Depending  on 
the  initial  conditions  either  one  of  these  equilibrium  positions  is  possible.  Fig.  O.b,  (3) 
Small  A,  and  moderate  velocity  parameter  A  will  give  simple  harmonic  response,  which 
correspond  to  elliptic  limit  cycle  in  a  phase  space.  Fig.  O.c,  (4)  The  more  complicated 
periodic  limit  cycle  motion  is  comprised  of  a  smaller  orbit  about  each  buckled  state 
and  a  larger  orbit  that  evolves  from  the  simple  periodic  motion.  Fig.  O.d,  (5)  Finally, 
a  chaotic  response  is  possible. 

It  is  interesting  to  note  that  we  observed  an  even  more  complicated  periodic  limit  cycle. 
For  the  value  of  the  in-plane  load  A«  =  -4.0ir’,  and  with  the  velocity  parameter  between 
A  =  119.275  and  A  =  120.496,  the  limit  cycle  of  the  form  presented  in  the  Fig.  10  has 
been  found.  Deformed  version  of  this  limit  cycle  has  been  observed  in  the  maneuvering 
case,  see  Fig.  20.n. 

Another  interesting  result  is  related  with  the  influence  of  the  initial  conditions.  Al¬ 
though,  our  attempt  was  not  to  study  the  initial-condition  problem,  basins  of  at¬ 
tractions,  or  fractal  basin  boundaries,  a  certain  section  of  initial  condition  space  was 
examined.  The  sequence  of  Figs.  ll.a  -  ll.e  shows  the  transient  from  the  chaotic  re¬ 
sponse  (Fig.  ll.a,  with  initial  condition  a\  =  0.084)  to  a  simple  limit  cycle  (Fig.  ll.b, 
oi  s  0.251),  and  then  through  a  period  3  motion  (Figs.  ll.c  and  ll.d,  ai  =  0.600,  to 
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1.600  back  to  chaos^  (Fig.  ll.e,  ai  =  2.000),  just  by  changing  the  initial  condition. 
Similar  behavior  have  been  observed  by  Moon  and  Li  (1985)  and  Pezeshki  and  Dowell 
(1987)  for  the  Duffing  equation. 

In  order  to  separate  the  influence  of  the  initial  conditions  from  the  maneuvering  the 
initial  conditions  are  taken  to  be  the  same,  oi  =  1.0  (corresponding  to  a  positive  velocity 
of  the  flrst  degree  of  freedom). 

4.2  Character  of  Solution  in  Presence  of  Maneuvering 

In  order  to  tmderstand  the  dynamic  problem,  it  is  helpful  to  think  in  physical  terms. 
The  buckled  plate  is  a  part  of  the  maneuvering  plane,  and  it  is  being  excited  by  the 
one  side  airflow.  Depending  on  the  value  of  the  load  factor,  and  other  parameters  one 
or  two  flxed  points  may  be  found.  The  calculations  show  that  three  possible  types  of 
stable  trajectories  may  occur.  There  are  the  two  limit-cycle  orbits  around  each  of  the 
static  equilibria,  and  one  global  limit-cycle  orbit  around  both  static  fixed  points.  Note 
finally,  that  there  is  also  the  possibility  of  chaotic  response.  Figure  12  shows  selected 
phase-plan?  plots  of  the  panel  response  as  well  as  theirs  power  spectrum  estimates.  For 
the  base  case  of  A  =  150,  and  Rg  =  -S.Oir^,  various  loads  factors  were  studied.  Figiire 
12. a  corresponds  to  a  non-maneuvering  case.  The  motion  has  non-chaotic  character 
with  three  closed  orbits  in  evidence.  A  load  factor  n  =  1.28,  transforms  the  response 
into  a  chaotic  state  (see  Fig.  12.b)  while  a  load  factor  n  =  3  suppresses  the  periodic 
character  of  the  motion,  the  response  is  a  fixed  point,  and  all  motion  ceases  (see  Fig. 
12.c).  The  most  remarkable  result  occurs  at  n  =  1.67,  when  the  maneuvering  has 
suppressed  one  of  the  three  orbits  that  appears  for  small  n  (see  Fig.  12.d). 

4.S  Deterministic  Behavior 

4.3.1  Convergence 

It  is  of  interest  to  compare  results  using  various  degrees  of  fireedom  (t.e.,  number  of 
modes).  Such  a  comparison  is  made  for  the  load  factor  n  =  2  in  Fig.  13  where  the  plate 
amplitude  of  the  limit  cycle  is  given  as  a  function  of  A  for  n/M  =  0.023,  P  =  0.0,  =  0.0 

for  2,  4,  6,  and  8  modes.  For  six  modes  of  freedom  the  solution  appears  converged, 
whereas  with  only  two  modes  the  results  are  inaccurate.  That  is  why  the  most  of 
the  calculations  presented  here  has  been  performed  using  six  degrees  of  fireedom.  Also 
shown  in  Fig.  13  are  non-maneuvering  results  for  2  and  4  modes.  Similar  behavior  of 
the  convergence  process  is  evident. 

4.3.2  Effect  of  Load  Factor 

In  Fig.  14  the  panel  amplitude  of  the  limit  cycle  is  given  as  a  function  of  A  for  several 
load  factors  n.  For  aeronautical  applications  n  is  typically  in  the  range  1 4-  6.  Several 

*  Word  of  csatino  ia  ^propriato  han:  the  roapon— a  ia  tha  Fisa.  ll.a  and  ll.a  mi^t  aa  waU  ba  quaaipariodic 
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trends  are  evident.  First,  the  flutter  speed  increases  with  n,  (i.e.,  the  linear  solution  is 
stabilized  by  the  presence  of  the  load  factor).  Secondlly,  the  amplitude  is  increasing  with 
increasing  the  load  factor.  As  might  be  expected  on  physical  grounds,  calculations  show 
an  increase  of  the  mean  amplitude  with  increasing  the  load  factor.  For  the  relatively 
large  n  the  response  is  a  fixed  point;  the  given  flow  velocity  is  not  large  enough  to 
disturb  the  static  equilibrium  shape  (see  also  Fig.  12.c).  Note  that  increasing  the  load 
factor  is  related  to  the  increase  of  the  static  pressure  P,  which  has  a  destabilizing  effect. 
Such  interaction  is  investigated  in  the  preceding  subsection. 

4.S.S  Effect  of  Static  Pressure  Differential 

Another  parameter  of  interest  is  the  uniform  static  pressure  differential  across  the  plate. 
In  the  absence  of  maneuvering,  the  plate  would  deform  to  some  static  equilibrium  po¬ 
sition  under  the  pressiure  differential.  With  maneuvering,  the  inertial  forces  due  to 
maneuvering  would  decrease  the  deformation  of  the  plate  due  to  static  pressure  differ¬ 
ential.  For  sufficiently  large  flow  velocity,  this  equilibrium  position  becomes  unstable. 
The  plate  then  begins  to  oscillate  about  an  unstable  equilibrium  shape,  eventually 
reaching  a  new  stable  dynamic  equilibrium  shape  (a  limit  cycle). 

The  results  of  the  analysis  of  the  influence  of  pressure  differential  to  the  amplitude  of 
panel  response,  for  a  given  load  factor  n  =  3,  are  presented  in  Figs.  15,  and  16.  In 
Fig.  15  the  amplitudes  of  the  periodic  response  are  plotted  as  function  of  dynamic 
pressure  A,  for  several  (nondimensional)  pressure  differentials  P.  It  u  apparent  that 
tl  e  flutter  speed  decreases  with  increasing  the  static  pressure  differential.  Furthermore, 
the  amplitudes  are  slightly  decreasing  with  the  increasing  of  the  static  pressure.  In  Fig. 
16  the  mean  value  of  the  amplitudes  of  the  limit  cycle  oscillation  are  plotted  as  function 
of  dynamic  pressure.  For  the  given  load  factor  n  =  3,  increasing  of  the  static  pressure 
differential  is  related  to  increasing  the  mean  amplitude  of  the  limit  cycle  oscillations. 

In  Fig.  17  the  influence  of  the  maneuvering  to  the  tnRyitniiin  and  minimum  values 
of  the  plate  deflection  is  given.  The  static  pressure  differential  was  constant,  P  =  60, 
during  thb  analysis.  The  presence  of  the  maneuvering  is  related  to  decreasing  of  the 
plate  deformation  caused  by  static  pressure  differential.  For  the  critical  value  of  the 
parameter  A  =  361  in  the  non-maneuvering  case,  and  A  =  353  in  the  maneuvering  case, 
a  supercritical  bifurcation  occurs.  The  fixed  point  in  the  parameter  space  bifurcate  and 
the  new  simple  harmonic  limit  cycle  motion,  with  lero  amplitude  at  the  bifurcation 
point,  is  created.  Note  that  in  the  range  A  =  353  -i-  361  the  presence  of  the  maneuvering 
changes  the  character  of  the  response  in  the  phase  space,  from  a  fixed  point  (n  =  l),  to 
a  limit  cycle  (n  =  2). 

4.S.4  Effect  of  In-Plane  Loading 

In  this  subsection  combined  infiuence  of  the  maneuvering  and  in-plane  loading  has  been 
studied.  Results  are  presented  in  Fig.  18.  The  plate  amplitude  is  shown  vs  dynamic 
pressure.  A,  for  several  load  factors,  n  =  1,2,  and  3,  and  in-plane  load  R,  =  -3/2ir’. 
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For  a  small  flow  velocities  A  there  are  two  branches  to  the  curve  associated  with  buck¬ 
led  configurations.  Depending  on  the  initial  conditions  either  one  of  these  nonlinear 
equilibrium  positions  is  possible.  With  maneuvering,  the  upper  buckled  position  ap¬ 
pears  only  for  very  small  A.  For  a  given  U,  =  -3/2ir*,  with  increasing  the  flow  velocity 
the  solution  is  stable  fixed  point  in  parameter  space.  Since  this  point  is  stable,  it  is 
an  attractor.  The  changing  in  the  nature  of  the  stable  solution  occurs  for  the  critical 
parameter  value  Ae  =  228  for  the  non-maneuvering  n  =  1,  and  A^  =  230  and  233  for 
the  maneuvering  conditions  n  =  2  and  3  respectively,  through  a  supercritical  or  normal 
bifurcation.  The  new  simple  harmonic  limit  cycle  motion  with  a  zero  amplitude  at  the 
bifurcation  point  is  created.  It  is  apparent  that  increasing  of  the  load  factor  is  related 
to  increasing  the  critical  parameter  value  A  at  the  bifurcation  point. 

4.4  Routes  to  Chaos 

In  the  following  we  will  discussed  how  the  maneuvering  transforms  the  response  from 
the  steady  state  regime  to  the  chaotic  regime.  The  attempt  is  to  answer  the  question 
how  is  a  chaotic  regime  established.  For  the  base  case  of  A  =  150,  and  A.  =  -3.0ir^,  load 
factor  was  varied  from  the  non-maneuvering  response  to  a  value  for  which  the  response 
is  a  fixed  point.  The  time  step  was  0.0025  (nondimensional  time),  the  number  of  time 
steps  were  60,000,  and  every  other  time  step  from  the  last  10,000  time  steps  has  been 
plotted.  Near  the  bifurcation  points  where  achieving  a  steady  state  solution  requires 
many  cycles  the  longer  time  interval,  90,000  time  steps,  has  been  used  to  avoid  the 
transient  regime.  In  order  to  decrease  possibility  of  mixing  the  numerical  instabilities 
with  the  chaotic  behavior  of  the  system  the  representative  cases  were  recalculated 
with  a  time  steps  of  0.001,  and  0.0005.  Figures  20.a  to  20.r  (discussed  in  details  in  the 
following  subsections)  shows  the  phase  portrmts  of  distinct  attractors  together  with  the 
power  spectrum  plots  discovered  during  the  numerical  experimentation.  Time  traces 
of  selected  responses  are  shown  in  Fig.  21. 

4.4.1  Intermittency  and  IVansient  Chaos 

We  begin  with  the  non-maneuvering  (n  =  l)  response  shown  in  Figs.  20.a,  and  21. a. 
The  motion  has  nonchaotic  character  with  three  closed  orbits  in  evidence.  For  the 
maneuvering  case  up  to  n  =  1.66,  we  see  what  appears  to  be  a  chaotic  response.  The 
power  spectrum  shows  large  fundamental  frequency  components  accompanied  by  a  less 
broad  band  noise  than  that  for  the  true  chaotic  response,  see  Figs.  20.b  and  20.c. 
Furthermore,  the  time  traces  show  the  periodic  motion  with  a  short  bursts  of  chaotic 
transient,  see  Fig.  21.c.  All  this  suggest  that  this  is  a  region  of  transient  chaos.  By  the 
value  of  the  load  factor  n  =  1.67,  a  new  stable  attractor  appears  (see  Fig.  20.d).  This 
triyectory  experiences  a  cascade  of  period-doubling  bifurcations  culminating  in  chaos. 
This  process  is  discussed  in  following  subsection. 
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4.4.2  Period-Doubling  Bifurcation 

A  period  doubling  bifurcation  is  easily  seen  for  n  =  1.83  (see  Fig.  20.e)  and  is  well 
established  for  n  =  1.90  (see  Fig.  20.f).  The  usual  splitting  of  the  trajectory  is  seen 
in  the  phase  portrait.  The  spectrum  also  shows  the  bifurcation  with  the  presence  of 
new  spectral  lines.  The  new  spectral  lines  show  that  the  period  has  now  doubled. 
Increasing  n  further  causes  the  previous  pattern  to  repeat;  a  splitting  of  each  of  the 
previous  harmonic  bands  into  two  bands  as  shown  in  Fig.  20.g.  One  anticipates  that 
by  continuing  to  increase  n,  the  same  phenomenon  will  be  repeated.  We  expect  to 
see  a  cascade  of  bifurcation,  each  accomipanied  by  the  period  doubling  associated  with 
a  Bubharmonic  instability,  leading  to  chaos.  The  numerical  calculations  support  this 
anticipation,  showing  the  chaotic  response  for  n  =  2.02. 

4.4.5  The  Inverse  Cascade  and  More  Chaos 

We  have  seen  that  the  approach  to  chaos  via  period  doubling  is  a  highly  structured 
process.  One  element  of  order,  to  be  observed  in  Figs.  20.i  thru  20.k,  is  the  existence 
of  normal  period  demultiplying  bifurcations,  culminating  in  another  new  limit  cycle, 
see  Fig.  20.k. 

Further  increase  in  n  cause  the  region,  n  =  2.44  -  2.97,  of  the  chaotic  response  to 
occur.  An  unexpected  limit  cycle  (at  n  =  2.66,  Fig.  20.n)  appears  in  an  island  bounded 
by  chaos  on  both  sides,  (Figs.  20.1  thru  20.q).  Uppon  further  examination  one  may 
cor.clude  that  this  is  basically  the  new  limit  cycle  discovered  in  the  non-maneuvering 
cate  (see  Fig.  10),  after  the  maneuvering  has  suppressed  one  of  the  orbits. 

Finally,  when  a  load  factor  reaches  n  =  3.00  the  inertial  forces  suppress  the  chaotic 
character  of  the  motion,  the  response  is  a  fixed  point,  and  all  motion  ceases,  see  Fig. 
20.r.  As  n  continues  to  increase,  the  defiection  of  the  plate  increases 

4.5  Nonautonomous  Equations 

In  the  case  when  the  angular  velocity  of  maneuvering  is  a  given  function  of  time,  the 
differential  equations,  Eq.  2.27,  are  nonautonomus  and  they  will  be  briefiy  considered 
in  this  subsection.  We  are  considering  an  airplane  in  a  puU-up  maneuver  during  the 
time  interval  -(3.  The  angular  velocity  is  given  by  w  =  mb  ■in(a<),  with  a  =  -<3). 

Under  these  assumptions,  the  angular  velocity  at  the  b^;inning  and  at  the  end  of  the 
maneuver  is  equal  to  sero  (see  Fig.  22).  As  might  be  e3q>ected  on  physical  grounds, 
our  calculations  show  a  disturbance  of  the  limit  cycle  by  the  inertial  forces  due  to 
maneuvering.  After  the  completion  of  the  maneuver  the  panel  returns  to  its  previous 
motion.  Phase  plane  response  of  the  panel  is  given  in  the  Fig.23. 

In  some  regions  of  the  initial  condition  space  one  may  e3q>ect  more  dramatic  changes 
of  the  dynamical  behavior  due  to  time  dependent  maneuver,  especially  close  to  the 
boundaries  of  the  basin  of  attraction.  It  is  clear  that  this  study  has  to  be  related  to  the 
analysis  of  the  infiuence  of  initial  conditions.  This  may  be  a  subject  of  future  research. 


SECTION  5 

CONCLUDING  REMARKS 

A  gene/al  geometrically-exact  Lagrangian  mechanics  formulation  for  the  aeroelastic 
analysis  of  a  maneuvering  aircraft  has  been  specialized  to  the  case  of  a  fluttering  plate 
imdergoing  a  pitching  maneuver.  The  formulation  includes  geometric  non-linearities 
associated  with  the  occurrence  of  tensile  stresses  in  the  middle  surface,  as  well  as 
the  effect  of  the  rigid-body  rotation  on  the  other  degrees  of  freedom.  The  general 
response  of  the  system  were  simulated  on  the  digital  computer  by  using  a  fourth- 
order  Rimge-Kutta  algorithm  with  adaptive  step  size  control.  Long-time  histories, 
three-dimensional  view  of  orbits,  phase  planes,  and  power  spectra  have  been  used  to 
characterize  the  response. 

We  have  chosen  system  parameters  such  that  there  are  two  equilibrium  position  for  the 
non-maneuvering  case.  The  numerical  simulation  for  periodic  and  chaotic  response  are 
conducted  in  order  to  analyze  the  influence  of  the  maneuvering  to  the  dynamic  behavior 
of  the  panel.  Several  important  conclusions  derive  from  the  preceding  work.  A  new  type 
of  limit  cycle  has  been  observed  in  the  non-maneuvering  case.  It  was  shown  that  chaos 
could  occur  in  a  maneuvering  case  for  system  parameters  in  the  actual  flight  range. 
The  presence  of  a  load  factor  n  can  transform  the  response  firom  the  flxed  point  into  a 
simple  periodic  or  even  chaotic  state.  It  can  also  suppress  the  periodic  character  of  the 
motion,  transforming  the  response  into  a  flxed  point,  and  all  motion  ceases.  Regions  of 
the  chaotic  motion  with  an  island  of  non-sinusoidal  periodic  motion  has  been  observed. 
When  n  =  1.83  we  observed  a  cascade  of  period-multiplying  bifurcations,  culminating 
in  chaos,  followed  by  a  series  of  period-demultiplying  bifurcations  as  the  system  returns 
to  a  starting  limit  cycle  response.  The  region  of  intermittent  and  transient  chaos  has 
been  observed  as  well. 

The  numerical  experiments  were  performed  with  a  different  time  simulation  length  in 
order  to  ensure  that  the  steady  state  response  has  been  reached  and  the  transient  has 
decayed.  Furthermore,  the  calculations  were  performed  with  various  step  size.  The 
results  were  in  good  agreement  and  showed  that  for  a  sufficiently  small  time  step, 
no  numerical  instability  occurred  and  the  results  for  the  time  simulation  were  closely 
reproducible. 

The  results  indicate  that  the  study  of  this  deterministic  system  is  important  from  the 
practical,  and  theoretical  viewpoint.  First,  the  fact  that  maneuvering  can  change  the 
character  of  the  panel  response  is  of  practical  interest  to  the  aeroelasticians  as  it  affects, 
for  instance,  fatigue  analysis.  Secondlly,  the  techniques  employed  in  this  study  can  be 
extended  to  the  problem  of  aeroelasticity  of  aircraft  (with  all  non-linearities  geometric, 
dynamic,  and  aerodynamica  included).  From  the  theoretical  point  of  view,  this  physical 


24 


system  is  a  rich  source  of  static  and  dynamic  instabilities  and  of  associated  limit-cycle 
motions  and  it  could  be  used  for  instance  as  a  test  case  for  assessing  techniques  for  the 
study  of  nonlinear  dynamics  and  chaos. 

Subjects  for  the  consideration  in  future  research  include  the  use  of  other  stability  con¬ 
cepts  such  as  Liapimov  exponents,  and  Poincare  maps.  Application  of  a  complementary 
methods  of  differentiable  dynamics,  in  particular,  of  center  manifold  and  bifurcation 
theory,  to  analyze  the  problem  from  a  qualitative  viewpoint  would  be  helpful.  Knowl¬ 
edge  of  the  generic  structures  of  attracting  sets  in  JV-space,  might  make  the  interpre¬ 
tation  of  numerical  solutions  of  evolution  equations  considerably  clearer. 
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Figure  1(a)  Airplane  in  a  puJl-up 
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Figure  1(b)  Budded  plale  of  a  mancuTcriiig  airplaae 
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Figure  3.  Diagrams  of  supercritical  bifurcations 

a)  Saddle-node 

b)  Transcritical 

c)  Pitchfork 

d)  Hopf 


S2 


Figure  4.  Diagrams  of  subcritical  bifurcations 

a)  Saddle-node 

b)  Transcritical 

c)  Pitchfork 

d)  Hopf 
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Figure  5.  Transition  from  a  stable  equilibrium 
to  a  strange  attractor 


Figure  6.  Collision  Figure  7.  Doubling 


Figure  8.  Birth  or  death  of  a  torus 
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Figure  9  Sketch  of  representative  phase  plane  orbits 
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PHASE  PLANE  PLOT 
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Figure  11(e)  Phase  plane  response  ai  =  0.084  (chaos) 
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Limit  cycle 
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Fignre  11(b)  Phase  plane  response  ai  =  0.251  (period  1  motion) 
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Figure  11(b)  C<mtiimed 


Figure  11(c)  Continued 
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FigoK  11(d)  Phaae  pUae  rcipoose  ai  =  1.600  (period  3  motioa) 
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Figure  11(d)  Contiiioed 
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PHASE  PLANE  PLOT 
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Figmc  12(d)  Contimied 
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CONVERGENCE  STUDY 


Figure  13  Amplitude  convergence  study 


MANEUVERING  EFFECTS 


Figure  14  Effect  of  load  foctor 
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Figure  IS  Ixifluence  of  static  pressure  differential 
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Figure  16  Influence  of  static  pressure  to  mean  amplitude  of  limit  cycle 
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PLANE  PLOT 


Figure  20  Continued 
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Figure  20  Continued 
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Figure  20  Continued 
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Figure  22  Airplane  in  time  dependent  poll-19  manenTer 
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